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Abstract 

Angular-energy distributions are studied for charged leptons and neutrinos from 
the decays of polarised top quarks. A small admixture of V+A interactions is 
incorporated. The polarisation dependent part of the neutrino distribution which 
can be measured experimentally through the missing momentum is particularly 
sensitive towards deviations from the V-A structure. This result remains unaf- 
fected by QCD corrections which, however, cannot be neglected in a quantitative 
analysis. 
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The analysis of polarised quarks and their decays, in particular of top 
quarks, has recently attracted considerable attention. Studies at a linear 
electron-positron collider are particularly clean for precision tests |T], |2| . How- 
ever, also hadronic [||, f|, U or 77 collisions || and subsequent spin analysis 
of top quarks might reveal new information. These studies will allow to 
determine the top quark coupling to the W and Z boson, confirming the 
predictions of the Standard Model or providing clues for physics beyond. 
For the top quark this latter possibility is particularly intriguing since it 
plays an exceptional role in the fermion mass spectrum. In this article we 
point out that neutrino momentum distributions are particularly sensitive 
towards violations of the V-A structure of the charged current. To perform 
such tests a sample of maximally polarised top quarks is required. 
A number of mechanisms have been suggested that will lead to polarised 
top quarks. For 77 collisions with circular polarised photons this possibility 
has been discussed in 0. Related studies may be performed in hadronic 
collisions which in this case, however, are based on the correlation between 
t and t decay products j|, 0, |5|. The most efficient and flexible reactions to 
produce polarised top quarks are electron-positron collisions. A small com- 
ponent of polarisation transverse to the production plane is induced by final 
state interactions which have been calculated in perturbative QCD ||, 
The longitudinal polarisation Pl is large. Its dependence on the production 
angle, beam energy and the top mass is discussed in [0, §. Pl varies strongly 
with the production angle, e.g. between nearly 0.6 for cos-# = —1 and —0.3 
for cos$ = 1 at y/s = 500 GeV. Averaging over the production angle leads 
therefore to a significant reduction of Pl with typical values of (Pl) around 
0.2 |J. QCD corrections change (Pl) by a relative amount of about 3% ||. 
All these reactions lead to sizable polarisation and can be used to obtain 
information on the production mechanism. The most efficient analyser of 
top polarisation is the charged lepton direction (as seen from the top rest- 
frame) in semileptonic top decays t — > Wb — ► l + vb. The decay distribution 
factorises into an energy and an angular dependent part |I| [ll] 

— (l + cos#+)/2 (1) 



dxi d cos 9 + dxi 

where 8 + denotes the angle between the top spin and the lepton direction. 
This factorised form applies for arbitrary top mass below and above the 
threshold for decays into real W bosons. QCD corrections to this distribution 



are well under control |T2[ and will be included in the discussion below. The 
analysing power of the charged lepton direction is obviously maximal and 
hence far superior to the W direction |I3J which suffers from the suppression 



factor (mj-2m^) / (mf +2m^). However, two drawbacks of all these proposals 
are evident: production and decay are mixed in an intricate manner, and 
furthermore the degree of polarisation is relatively small and depends on 
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the production angle. Top quark production with longitudinally polarised 
electron beams and close to threshold provides one important exception: the 
restricted phase space leads to an amplitude which is dominantly S-wave 
such that the electron (and positron) spin is directly transferred to the top 
quark. Close to threshold and with longitudinally polarised electrons one 
may deal with a highly polarised sample of top quarks independent of the 
production dynamics. Thus one may study the V-A structure of t decays 
under particularly convenient conditions: large event rates, well identified 
restframe of the top quark, and large degree of polarisation. 
As stated before the angular distribution of charged leptons is optimal for the 
analysis of the top polarisation. On the other hand exactly for this reason it 
is less suited to identify small admixtures of non-standard couplings, e.g. a 
small V+A amplitude. In the analysis of neutrino (= missing momentum) 
distributions the situation is reversed: neutrino distributions in top decays 
are sensitive to small V+A admixtures and less sensitive to top polarisation. 
This observation and the detailed analysis of charged lepton and neutrino 
distributions with a small admixture of V+A interaction is the key point 
of our paper. Two aspects will be considered. The strong sensitivity of 
the neutrino energy- angular distribution towards a V+A admixture will be 
demonstrated. In addition it will be shown that QCD corrections to the 
Standard Model distributions are well understood, such that any deviation 
can truly be attributed to a deviation from pure V-A coupling in the decay. 

As stated before, the emphasis of this work will be on a discussion of po- 
larised top quarks decays in their rest frame. Nevertheless, a brief comment 
is appropriate on relativistic top quarks produced far above threshold at a fu- 
ture linear collider. With increasing energy it will be more and more difficult 
to reconstruct the double differential energy- angular distributions for the lep- 
tons. For ultra-relativistic top quarks the situation will resemble the one for 
b quark semiloptonic decays at LEP energies. A difference which makes the 
studies at the Z° peak particularly interesting is the very high polarisation 
of b quarks, Pl = (—0.94), which weakly depends on the production angle 
[H]. For top quarks the net polarisation will be smaller, as stated earlier, 
but evidently a lot may be gained when considering the asymmetry or the 
first Legendre moment rather than the average over the production angle. 
Moreover, the short lifetime in the top case helps to avoid depolarisation in 
hadronization which is a serious problem for the bottom quark. Thus the 
studies for ultra-relativistic top quarks are also worth consideration. In close 
analogy to b quarks at LEP one can study the charged lepton energy distri- 



bution (and its moments) in the laboratory frame. Formulae are given in 



which relate this distribution to the double differential angular-energy dis- 
tribution discussed in the present article. The neutrino distribution, which 
can be treated in the same way, has recently attracted considerable inter- 
est. It has been pointed out |l5j that for b decays the neutrino distribution 
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is highly sensitive to the polarisation of the parent quark and can be used 
in polarisation studies for A& baryons in addition to the distribution of the 
charged lepton. In it has been advocated that the ratio of the average 



energies of the charged lepton and the neutrino is particularly sensitive to 



b polarisation. In a recent preprint (17[ a study has been proposed of the 



charged weak current and its space-time structure in b semileptonic decays 
through the neutrino distribution. Unlike for the top decay where the in- 
variant mass of the lepton system is equal to m w , the mass of W boson, 
for b decays one usually integrates over this variable. We think, however, 
that another difference is much more important. Neglecting the effects of W 
propagator the neutrino spectrum in the decays of a heavy quark with the 
weak isospin ^3 = —1/2 is the same as the spectrum of the charged lepton for 
^3 = 1/2. Thus the results of the present article suggest that for polarised 
b quarks the distribution of the charged lepton should be more sensitive to 
V+A admixture. The detailed discussion will be given elsewhere. 

Let us perform now a quantitative discussion of the sensitivity of different 
distributions to V+A admixtures. The tbW vertex is parametrised as follows: 

with 



g v = (l + «)/vT+^, 9a =(-1 + k)/VTT^ (2) 

Hence k = corresponds to pure V-A and k = 00 to V+A. We consider 
only leptonic decays of W. The masses of all fermions in the final state are 
neglected. The error associated with the lepton masses is negligible at least 
for electrons and muons. The effects of non-vanishing b quark mass are of 
order ne, where e = mf,/m t , so they are ~ 1% for k = 0.3 and e = 0.03. The 
following discussion refers to the rest frame of the decaying top quark. 
Neglecting QCD corrections the differential decay rate is proportional to the 
following expression: 

dr l — (b- vR_ ■£ + n 2 b-lR+ ■ u) (3) 

1 + K ^ J W 

where 

R± = t±s 

t, b, i and v are the scaled four-momenta of top, bottom, charged lepton and 
neutrino. The common scaling factor is mj~ , so 

t 2 = 1 and y = (£ + z/) 2 = (m w /m t ) 2 . 

s — (0, s) is the spin four- vector of the decaying quark. 5* = \s\ = 1 corre- 
sponds to fully polarised and S = to unpolarised top quarks. The following 
variables will be used: 



xi = It- 1 Scos6 + = -s-l/t-l 

x v = It ■ v S+os^o = —s-vjt-v 



(4) 
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Thus, %t and proportional to the energies of the respective leptons, 9 + 

is the angle between s and the three- momentum of the charged lepton, and 6q 
is the analogous angle for the neutrino. Employing the Dalitz parametriza- 
tion of the phase space and following the steps described in Appendix B of 
|y| and Appendix A. 2 of [|Tj| one derives the following formulae for the nor- 
malized double differential angular and energy distributions for the charged 
lepton and for the neutrino 



diV 



dx/dcos^j 
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1 + ^) My) 
.2 



S cosd + J^(x £ ,y) 



r-r{F+(x g ,y) 

o(x t ,y) - 5 cos 6> + J (x e ,y)] } 

(y<x e <l) (5) 



dN 



dx u d cos 9q 



where 
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;i + My) 



{ F (x v , y) + S cos 6 J (x v , y) 



{x v ,y) - ScosOq Jo (x v ,y) } 



K&y) 
Jo ( x >y) 
' 



Fn (x,y) 



J ( x >y) 
My) 



x(l — x) 

K&v) 

(x-y)(l -x + y) 

(x - y)(l - x + y - 2y/x) 

2(l-y) 2 (l + 2y) 



(y<x u <l) (6) 



(7) 
(8) 
(9) 
(10) 

(11) 



These distributions can be easily integrated over the energies leading to the 
following angular distributions 



dN 



d cos 9 + 

dN 
d cos 6Y1 



— I — S cos 8 j 
2 2 H 



1 



1 + k' 



Ky) 



1 
2 



5 cos O [1 - %)] 



1 + 



K 



Ky) 



1 + K 2 1 - %) 



where 



%) 



123/(1 - y + 2/lnj/) 



(12) 
(13) 

(14) 



(1 - 2 /) 2 (l + 2 2/ ) 

For m w = 80 GeV and m t = 160 GeV one obtains y = 0.25 and h(0.25) = 
0.566. It is evident that the neutrino angular distribution is significantly 
more sensitive towards the admixture of V+A than the angular distribution 
of the charged lepton. One can also gain sensitivity by studying the double 
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differential angular-energy distributions. A convenient method is to calculate 
moments of the distributions. The average energies provide the simplest 
example: 



d(x £ ) 
d cos 9-> 



d(x u ) 
d cos 9 



l + 2y + 3y 2 
4(1 + 22/) 

+ S cos 9 + 

l + ^y + y 2 
4(1 + 2y) 

i 

+ S cos 9 



1 + 



^ 2y(l - y) 
1 + k 2 l + 2y + 3y 2 
k 2 2(1 -3y + 2y 2 ) 



1 + k 2 l + 2y + 3y 2 
k 2 2y(l-y) 



(15) 



8y + y 



1 + /«' 

2 



4(1 + 2^/) 



2 1 
1 - 



4y + 2/ 2 
k 2 2(1 



3y + 2y 2 



K 



l-8y + y 2 



(16) 



Hence the V-A form of the charged current in top decays can be tested by 
measuring these averages. Even without polarisation they are changed by 
an admixture of V+A, so, in principle at least, it is possible to see some 
effects although the 'signal/background' ratio is smaller. The best strategy, 
however, is to combine the results on both the charged lepton and the neu- 
trino distributions for highly polarised top quarks. At this point it becomes 
necessary to include QCD corrections to the decay. 

QCD corrections to the lepton spectra in decays of polarised quarks were 
calculated previously [] |y], [12|, [L5|, [T5| . It is plausible that they reduce the 



V+A contribution to the total rate by the same factor as the well-known 

p0|1 . This implies that the ratio of V-A and V+A 
also after inclusion of 0(a s ) corrections. This 



one for V-A contribution 
contributions should be 1 : k v 
condition will be imposed on the QCD corrected normalized distributions 
given in the following. Anyway, the unknown corrections are small (~ a s n 2 ) 
and the errors related to our procedure should be even smaller. The normal- 
ized distributions including first order QCD corrections read: 



diV 



dxfdcos9. 



dN 



dx u d cos 6*o 



1 
2 
1 
2 



A e (x e ) + Scos9 + B e (xe) 
,&- v {x v ) + S cos 9 B 1/ (x l/ ) 



(17) 
(18) 



lr The results of are in conflict with earlier analytic calculations . The QCD 
correction to the lifetime of t quark |2(J follows from and has been confirmed by other 
groups pH ■ The result of pd| ] on the e + spectrum in charm decays has been also confirmed 
p^ |. For e~ from b decays see Q, footnote on page 27. Adding this to the cross-checks 
performed in [JLlJ [1^] the controversy can be considered as solved in favor of our article 
PH and the related subsequent articles. 
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Ad 


r x) = 


12 


~K(x,y) 


- a s Ff(x,y) 


+ 


k 2 Fo(x,y) 


1 + K 2 


My) 


- a s My) 


My) 


Be 




12 




- a s Jt(x,y) 




k 2 Jo(x,y) 


1 + K 2 


My) 


- a s Mv) 




My) 


A 




12 


Fo{x,y) 


- a s F{(x,y) 






1 + 


My) 


- a s Mv) 


My) 


B u 


[x) = 


12 




- a s J{(x,y) 




« 2 Jt{x,y) 


1 + K? 


{ My) 


- a s Mv) 




My) 



(19) 
(20) 
(21) 
(22) 



where 



2a. 
3tt 



My) = My)[^ 2 + 4U 2 (y) + 2]nyln{l-y)]-{l-y)(5 + 9y 

+ 4y(l -y- 2y 2 ) \ny + 2(1 - y) 2 (5 + Ay) ln(l - y) 
as first derived in [0, and [JIT], 0, |TJ 

F+(x,y) = F+(x,y)$ + z$i-(3 + 2x + y)$ 2o3 + 5(l-x)$ 4 + (- 
+9x - 4a; 2 - 2y - y 2 )$ 5 + y(4 - 4x - y + y/z) /2 



-6y 2 ) 
(23) 

-2xy 
(24) 



Fi(x j3 /) = Fo(x,y)$ + (-2xy + x + y + y 2 )$ 1 + (-5 + 2a;-3y)$ 2<> 3 
+ (5 + 4xy -5x + 3y - 5y 2 - 2y 2 /x)<§> 4 + (Gxy + 9x - 4a; 2 
-lly - 2y 2 + 2y 2 /a;)$ 5 + y{2 + 3x - 3y - 2y/x)/2 (25) 



J^(x,y) = J f (x,y)$ -x$i + (5 - 2x - ?/ - 2y/x - 2/x)$ 2 o3 + (-3 + x 
+2/x)$ 4 + (-2zy + 3x - 4x 2 + 6y - y 2 - 2y 2 /x)§ 5 
+ (2 - 2a; 2 + 2y - 3y 2 - 2y/a; + 3y 2 /x)/2 (26) 



J7/(z,y) = ^{x,y)$Q + (-2xy-x-hy + y 2 - 2y 2 /x)$i + (3 + lOx + ?/ 
+10y/s - 2/a;)$ 2 o3 + (-3 + I2xy + x - 7y - y 2 - 12y/x 
+8y 2 /x + 2/x)$ 4 + (Gxy - 9x - 4x 2 - y - 2y 2 + 10y 2 /x)$5 
+ (2 - 5xy - 2x 2 + 2y + 7y 2 -2y/x- 2y 2 /x) /2 (27) 

where 

$o = ^ + 2Li 2 (x) + 2Li 2 (y/x) + ln 2 (i^) 
$i = ^ + Li 2 (y) -Li 2 (x) -Li 2 (y/x) 
$203 = |(l-j/)ln(l-y) 
$ 4 = |ln(l— x) 

$ 5 = |ln(l-y/z) (28) 
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It is worth mentioning that in [|T^] formulae are also given for the QCD 
corrections to V-A heavy quark decays for non-zero mass of the quark pro- 
duced in the decay which are only slightly more complicated than those in 
the massless approximation. Thus it is easy to generalize the present dis- 
cussion to non-zero e. As stated before, terms ~ ne appear for e ^ 0. The 
corresponding QCD corrections of 0(a s ne) are unknown but small, and can 
be safely neglected. Anyhow, in the present article we stick to the massless 
approximation for the sake of simplicity of the presentation. 



k 


4" 


Bf 


aV 




-1 


2.008 


2.005 


1.593 


-0.707 




.981 


.940 


1.023 


1.166 





1.000 


.998 


1.000 


-0.452 




1.000 


.949 


1.000 


1.110 


1 


.559 


.558 


.683 


-0.322 




1.021 


.960 


.984 


1.068 


2 


.345 


.344 


.500 


-0.249 




1.043 


.973 


.973 


1.037 


3 


.230 


.230 


.385 


-0.203 




1.064 


.989 


.965 


1.015 



Table 1: Moments of the coefficient functions defining the angular-energy 
distributions for the charged leptons and for the neutrinos for y = 0.25 and 
a s = 0.11: the upper entries denote the moments for k 2 = and the lower 
ones are the ratios between the moments for k 2 = 0.1 and n 2 = 0. 

In Table 1 the following moments are given 

Af = Jy dy x k A e (x, y) , Bf = ft dy x k B e (x, y) (29) 

A[ u) = fZdyx k A v (x,y) , = J y ' dy x k B u (x,y) (30) 

for integer k between -1 and 3, y = 0.25 and a s (m t ) = 0.11 . The upper 
entries in the table denote the values of the moments for n 2 = and the 
lower ones the ratios of the moments evaluated for k 2 = 0.1 to those for 
Kj 2 = 0. From the comparison of the moments it is again evident that the 
moments B% which govern the angular dependence of the neutrino spectrum 
are particularly sensitive towards a V+A admixture. The effect is most 
pronounced for the moment k = — 1 which enhaces the lower energy part of 
the spectrum and where the relative change amounts to 17% for k 2 = 0.1 . 
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Figure 1: The coefficient functions a) A^(x) and b) Bi(x) defining the charged 
lepton angular-energy distribution for y = 0.25 and a s {m t ) = 0.11 : k 2 = 
- solid lines and k 2 = 0.1 - dashed lines. 




Figure 2: The coefficient functions a) A u (x) and b) B u (x) defining the neu- 
trino angular-energy distribution for y = 0.25 and a s {m t ) = 0.11 : k 2 = 
- solid lines and k 2 = 0.1 - dashed lines. 
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The same conclusions follow from Figs.|T] and where the coefficient functions 
defined in eqs. (|l9l) -(|22|) are shown as solid lines for k = and as dashed lines 
for ac = 0.1 . In Figs.^a and 0a the functions Ag(x) and A u (x) are plotted 
for y = 0.25, a s (m t ) = 0.11 and the functions B^(x) and B„(x) are shown in 
Figs.|l]b and ||b, respectively. 





K 2 


A,(0.5) 


B,(0.5) 


A„(0.7) 


B„(0.7) 


0.11 





1.816 


1.812 


1.794 


-0.583 


0.11 


0.1 


1.772 


1.688 


1.767 


-0.666 








1.778 


1.778 


1.760 


-0.526 





0.1 


1.737 


1.657 


1.736 


-0.614 



Table 2: Comparison of the coefficient functions with and without QCD 
corrections and V+A admixture for y = 0.25 

The effect of QCD corrections is illustrated by the comparison of the coeffi- 
cient functions A^(0.5), B^(0.5), A^(0.7) and B u (0.7) which are given in Table 
2. The effect of QCD correction can mimic a small admicture of V+A in- 
teraction. Therefore, inclusion of the radiative QCD correction to the decay 
distributions is necessary for a quantitative study. 

We conclude that the angular-energy distribution of neutrinos from the 
polarised top quark decay will allow for a particularly sensitive test of the 
V-A structure of the charged current. 
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